Astronomy <fe Astrophysics manuscript no. Hgl82 


February 5, 2008 


(DOI: will be inserted by hand later) 





Spitzer IRS spectroscopy of 3CR radio galaxies and quasars: 

Testing the unified schemes 

Martin Haas 1 , Ralf Siebenmorgen 2 , Bernhard Schulz 3 , Endrik Kriigel 4 , and Rolf Chini 1 



in 
o 
o 

(N 

Or 

00 

m 



> 
o 

m 

On 
O 

in 
o 

Or 

o ■ 
a: 



1 Astronomisches Institut, Ruhr-Universitat, Universitatsstr. 150 /NA7, D-44780 Bochum, Germany 

2 European Southern Observatory, Karl-Schwarzschildstr. 2, D-85748 Garching b. Miinchen, Germany 

3 IPAC, California Institute of Technology (Caltech), Pasadena, CA 91125, USA 

4 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, Postfach 2024, D-53010 Bonn, Germany 
xx.07.2005/xxx 

Abstract. With the Spitzer Space Telescope Infrared Spectrograph (IRS) we have observed seven powerful FR2 
radio galaxies and seven quasars. Both samples, the galaxies and the quasars, are comparable in isotropic 178 



MHz luminosity (10 2b5 W/Hz < Pi 



< 10 29 E> W/Hz) and in redshift range (0.05 < z < 1.5). We find for both 



samples similar distributions in the luminosity ratios of high- to low-excitation lines ([NeV]i4.3 Mm / [NeII]i2.8pm) 
and of high-excitation line to radio power ([NeV]24.3fjm / Pi78 MHz). This solves the long debate about the apparent 
difference of quasars and radio galaxies in favor of the orientation-dependent unified schemes. Furthermore, the 
luminosity ratio [OIII]5oo.7nm / [OIV]25.9/im of most galaxies is by a factor of ten lower than that of the quasars. 
This suggests that the optical emission from the central NLR is essentially absorbed (Av ^ 3) in the powerful FR 2 
galaxies and that the [OIII]5oo.7nm luminosity does not serve as isotropic tracer for testing the unified schemes. 

Key words. Galaxies: active - Galaxies: nuclei - Galaxies: quasars: general - Infrared: galaxies 



1. Introduction 

According to the orientation-dependent AGN unified 
schemes both the powerful FR2 radio galaxies and the 
quasars are believed to belong to the same parent popula- 
tion (Barthel 1989) with their observed properties depend- 
ing on our viewing angle towards the obscuring dust torus 
(Antonucci 1993). First evidence for this hypothesis was 
provided by spectro-polarimetric observations in those 
galaxies with "edge-on" torus, where scattering particles 
located above and below the dust torus allow for viewing 
the AGN-typical broad spectral lines of the central region 
(e.g. Antonucci 1984). In general a proper test of the unifi- 
cation is best performed by observations of isotropic quan- 
tities, which are not affected by orientation-dependent ex- 
tinction or beaming. Since the fraction of the nuclear lumi- 
nosity that is absorbed by dust must be re-radiated in the 
infrared, the optically thin far-infrared (FIR, 40-1000 /im) 
emission is such an isotropic measure. From sensitive pho- 
tometry with the ISO satellite one found similar isotropic 
FIR dust luminosities for radio galaxies and quasars from 
the 3CR catalogue, after normalisation by their likewise 
isotropic 178 MHz radio power (Meisenheimer et al. 2001, 
Haas et al. 2004). The distributions showed a large dis- 
persion, which may be caused by the various evolutionary 
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stages of these complex sources, and it is still under debate 
whether the FIR luminosity arises from starbursts rather 
than from the AGN itself. 

Both the radio galaxies and the quasars are luminous 
mid-infrared (MIR, 10-40 (im) emitters, and the mid- to 
far-infrared luminosity ratio of the 3CR sources exceeds 
that of pure starburst galaxies; but the radio galaxies 
show a lower Lmir / Lfir than quasars (Haas et al. 2004, 
Siebenmorgen et al. 2004). Either the emission from the 
putative dust torus is optically thick even at MIR wave- 
lengths (Pier & Krolik 1993) - hence not isotropic -, or the 
intrinsic AGN luminosities are different supporting mod- 
ifications of the orientation-dependent unification such as 
the receding torus model (Lawrence 1991). Further objec- 
tions against the pure orientation-dependent unification 
come from the lower [OIII]5oo.7nm / Pi78MHz luminosity 
ratios of galaxies, assuming that most of the [OIII]5oo.7nm 
emission arises from the extended narrow line region not 
affected by central extinction. 

The high- and low-excitation MIR spectral lines suffer 
from much less extinction than optical ones. Therefore, 
we observed 7 quasars (and BLRGs) and 7 radio galaxies 
from the 3CR catalogue. First results from this sample on 
the discovery of 10 /im silicate emission in quasars have al- 
ready been presented by Siebenmorgen et al. (2005). Here, 
we focus on the investigation of the emission lines. 
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Fig. 1. IRS Spectra of the radio galaxies (top panels) and quasars (bottom panels), shifted to the object rest frames. 
The important emission lines as well as the silicate absorption and emission features are marked. For 3C047 and 
3C 445 only two of the four IRS channels were assigned to our program. 
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Fig. 2. Luminosity ratios of high over low excitation 
lines versus mid-to far-infrared luminosity. For 3C 047, 
3C079, 3C321 and 3C445) [Nell] was substituted by 
[S III], and for 3C047 [NeV] 14 . 3tim by [NeV] 24 .3 M m- The 
cross indicates errors of a factor 1.5 (MIR/FIR) and 2 
([NeV]/[NeII]). 



Fig. 3. High-excitation line to radio luminosities ver- 
sus [OIII] 50 o.7nm / [OIV] 2 5.9 M m luminosities. For 3C295, 
3C298 and 3C356 [NeV] 2 4.3 Mm was substituted by 
[NeV]i4. 3 pm, and for 3C047, 3C295, 3C298, 3C351 and 
3C356 [OIV] 25 .9p m by 3 x [NeV]. The cross indicates er- 
rors of a factor 2. 



2. Sample selection and Observations 

The steep radio spectrum quasars and powerful galaxies 
were randomly selected from the 3CR catalogue, so that 
the samples match in 178 MHz luminosity as well as in red- 
shift. Note that we do not consider FR 1 or low-power FR 2 
sources, since for them the orientation-dependent unifica- 
tion picture may not apply. Thus, any conclusions drawn 
here refer to FR2 sources with 10 26 5 W/Hz < Pi 78M Hz 
< lO 29 5 W/Hz and 0.05 < z < 1.5. 

More specifically, the sources were selected from that 
sub-sample of the 3CR catalogue which (1) was not re- 
served by guarantee time observations and (2) was ob- 
served by ISO with sufficient S/N (Haas et al. 2004, 
Siebenmorgen et al. 2004) in order to provide the flux 
estimates required for carefully planning the Spitzer IRS 
observations. This "ISO subset" provides a fairly represen- 
tative sub-sample of the 3CR catalogue, since the objects 
varied widely in apparent brightness and the selection cri- 
teria for the sources were low cirrus foreground and good 
visibility to the satellite. Thus, we do not expect that the 
targets selected for our Spitzer IRS observations are biased 
in favour of MIR-loud 3CR sources. 

The objects were observed between 5 and 35 /^m with 
the infrared spectrograph (IRS, Houck et al. 2004) of the 
Spitzer Space Telescope (Werner et al. 2004) in the two 
IRS low-resolution (64 < A / AA < 128) modules in star- 
ing mode. Our analysis starts with the two dimensional 
data frames from the Spitzer pipeline (Higdon et al. 2004) 
using the latest calibration files. At this point, the major 
instrumental effects have already been removed. We sub- 



tracted the sky background, using pairs of frames, where 
the source appears at two different positions along the 
spectrometer slit. We interactively extracted the one di- 
mensional spectra. Before averaging, parts with low repro- 
ducibility between integrations were discarded. 

3. Results and Discussion 

Figure ^] depicts the IRS spectra shifted to the object rest 
frames. The spectra show numerous emission lines as su- 
perposed on diverse continuum shapes, caused by various 
hot and warm dust components, and the broad silicate 
features around 10 and 18 /jm. The silicate features are 
typically seen in absorption for the galaxies (e.g. in 3C 079, 
3C 303.1, 3C321, 3C459) and in emission for the quasars 
(e.g. in 3C109, 3C 249.1, 3C 323.1, 3C351). Also in the 
NIR (3-10 fim) shortward of 7 /im the continuum shows a 
typical difference between galaxies and quasars: in quasars 
it is steeply rising with decreasing wavelength while in 
galaxies it is flat or declining. From both the different sil- 
icate features and the bluer NIR continuum of quasars we 
suggest that the reported difference of the Lmir/Lfir ra- 
tio (Lio-40/j m /L4o-iooo/xm) is in fact caused by suppression 
of the MIR continuum in galaxies and by its enhancement 
in quasars. 

The spectra show AGN-typical high-excitation lines 
like [NeV]i4.3^m, [NeV] 24 .3 M m and [OIV] 25 . 9pm , as well 
as low-excitation lines like [NeII]i 2 .8/j m and [S III] is. 7pm- 
Tabled lists the line luminosities. We estimate the uncer- 
tainty of the line fluxes from the subtraction of the un- 
derlying noisy continuum to ^40%; combined with the 
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Table 1. Sample parameters. The types are abbreviated as Q=quasar, B = broad line radio galaxy, N = narrow line 
radio galaxy. The luminosity distance Dl is computed adopting a A-cosmology with H = 70 km s -1 Mpc _1 , Q m = 0.27 
and fl\ = 0.73. Below the line wavelength we list the ionisation potential needed to create the species from the next 
lower ionisation stage. The line luminosities are given logarithmically in units of [W]. The [OIII]5oo.7 n m luminosities 
are taken from the compilation by Grimes et al. (2004) and from Boroson & Oke (1984), Tadhunter et al. (1993), Mc 
Carthy et al. (1995) and Hirst et al. (2003). 
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a line between two channels, b line on edge of channel, 



c adopting L([OIII]) — 0.1 X L(total H a ) from the mean SDSS quasar template (Vanden Berk et al. 2005) 



absolute flux calibration uncertainty of ^ 20% this results 
in a total line flux uncertainty of ~ 50% which is sufficient 
for the conclusions drawn here. In order to compare the 
two samples by means of representative lines, when such 
a line is unobserved or poorly measured, we substitute it 
by another one of similar (or higher) ionisation potential 
with appropriately scaled strength determined from aver- 
age line ratios: [NeII]i 2 .8 M m by [S III]i 8 . 7Aim , [NeV] 2 4.3Atm 
by [NeV]i4. 3 A«n (and vice versa), and [OIV] 2 5.9Ami by 
3 x [Ne V]24.3A 1 m- This leaves the scatter in the line ratios 
below a factor of two. 

Figure |U displays the [NeV] 14.3^1 / [Nell] 12.8pm lumi- 
nosity ratio versus Lmir/Lfir- While the radio galax- 
ies show on average a lower Lmir / Lfir value than 
the quasars (by a factor ^3), the distribution of 
[NeV] i4.3aoti / [Nell] i2.8Aim is strikingly similar for both 
samples: The logarithmically determined mean of the 
[NeV] / [Nell] ratios (and in bracketts the rms expressed 
as factor around this mean) are 1.1 (2.8) and 0.9 (2.1) 
for the quasars and galaxies, respectively. Furthermore, 
the [NeV] i4.3Aim / [NeII]i2.8/im values are much higher (by 
a factor ~100) than those found for starburst galax- 
ies (Sturm et al. 2002), placing also the FR2 galaxies 
clearly in the AGN dominated range. We find also similar 
results using other lines, e.g. [NeVI] 7 . 6/ml , [NeIII]i 5 . 5 ^ m , 
[SIV] 

10.5Atm ■ 

If the unification between powerful galaxies and 
quasars is correct, then the ratio of central AGN lumi- 
nosity to isotropic radio lobe power should be similar for 
both types. Fig. |31 (vertical axis) shows the luminosity ra- 
tio [NeV]24.3Aim / Pi78MHz- Again, the distributions look 
similar for both samples: The logarithmically determined 
mean of the ratios (and in bracketts the rms expressed as 
factor around this mean) are 0.027 (2.2) and 0.023 (4.1) 
for the quasars and galaxies, respectively. This suggests 
that the galaxies possess intrinsic AGN luminosities simi- 



lar to those of the quasars. It provides direct spectroscopic 
evidence in favour of the unified schemes. 

How can these results be reconciled with the fact found 
by many authors, that galaxies show a ten times lower 
[OIII] 500.7nm /Pi78MHz luminosity ratio than quasars? 
Does optical extinction play a role as proposed by e.g. 
Jackson & Browne (1990), Hes et al. (1996), Baker (1997) 
and di Serego Alighieri et al. (1997), or do galaxies possess 
a weaker central powerhouse suggested by e.g. Lawrence 
(1991) and Grimes et al. (2004)? We can now test these 
two possibilities using the [OIV]25.9Atm line: it has a higher 
ionisation potential than the optical [OIII] 500.711m line, but 
is ~50 times less affected by extinction. If the galaxies 
had intrinsically weaker AGN and hence weaker radiation 
fields, then one expects their [OIII] 5 oo.7nm / [OIV]25.9Aim 
ratios to be higher than for quasars. However, this is 
not the case, rather we find the opposite as shown in 
Fig. 13 (horizontal axis). The galaxies exhibit a lower 
[OIII]5oo.7nm / [OIV]25.9Aim luminosity ratio than all but 
one of the quasars. For this one quasar (3C 298) the 
[OIII]5oo.7nm is highly uncertain, since extrapolated from 
the total (broad and narrow line) H Q luminosity (Tab.0. 
The logarithmically determined mean of the [OIII] / [OIV] 
ratios (and in bracketts the rms expressed as factor around 
this mean) are 16.9 (2.2) and 2.5 (2.2) for the quasars and 
galaxies, respectively, hence the distributions are different 
at more than 5c Therefore we suggest that the galaxies 
suffer from substantial [OIII]5oo.7nm extinction (Ay <^3), 
at least in the central regions and that most of the - in 
principle extended - NLR [OIII]5oo.7nm luminosity origi- 
nates from the central part. This result questions the for- 
mer use of the [OIII]5oo.7nm emission as isotropic tracer 
for testing the unified schemes, in particular for establish- 
ing the receding torus model. As proposed by Haas et al. 
(2004), it seems that not only the inner radius, but also 
the height of the torus increases with growing luminosity. 
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In contrast to our high luminosity 3CR objects, for 
Scyfert galaxies the luminosity ratio [OIII] 5 oo.7nm/FIR is 
the same for type 1 and type 2 when the objects are se- 
lected by LpiR with a 25/im/60^im colour criterion (Keel 
et al. 1994, Fig 3). This suggests that the [OIII] 50 o.7nm 
emission is isotropic for those objects of relatively low 
luminosity. On the other hand, the mostly radio-quiet 
Sloan AGN seem to indicate that the ratio of type Is 
to type 2s increases with luminosity, when selected by 
[OIII]5oo.7nm (Simpson 2005). This result can be qualita- 
tively explained in the context of the unified model with- 
out recourse to a luminosity-dependent torus opening an- 
gle if the [OIII]5oo.7nm emission is anisotropic at high in- 
trinsic optical/UV luminosity. For the radio galaxies, if 
the [OIII] 5 oo.7nm emission is mostly hidden, then relative 
to the broad lines it should be almost as strong in polarised 
flux as in total flux. Notably in three of five narrow-line 
radio galaxies (0.11 < z < 0.22), di Serego Alighieri et al. 
(1997) found polarised [OIII]5oo.7nm emission. 

One source (the anomalous object 3C459) has only 
rather high upper limits in the high-excitation [NeV] lines, 
which are still consistent with the unification. Its contin- 
uum is probably mostly starlight from a post-starburst 
population (Miller 1981). Also, though classified as an 
NLRG here, some authors claim to have detected broad 
lines (see the "notes" in NED). Notably, its IRS spec- 
trum exhibits excited molecular hydrogen emission (H 2 
S(l) at 17.03/j.m and H2 S(3) at 9.67/zm with luminosi- 
ties of 10 33 - 78 W and 10 3434 W, respectively) - to our 
knowledge the first such detection in a powerful FR 2 ra- 
dio galaxy. Compared with [OIV] 25.9^111 and Pi78MHz this 
source shows also an FIR luminosity excess (Haas et al. 
2004). Since 3C459 is a compact steep spectrum radio 
source, we suggest that the excited molecular hydrogen 
emission originates from shocks due to the probably young 
expanding jet. 

On the other hand, it is worth to note that the 
radio galaxy 3C295 with a very low Lfir/Pi7smhz 
(Meisenheimer et al. 2001) lies still at the low end of the 
[Ne V] / Pi78MHz distribution, but has a rather interme- 
diate [NeV] /[Nell] ratio; this indicates that despite the 
AGN-typical line ratio 3C 295 has exceptionally bright ra- 
dio lobes possibly caused by a special evolution or envi- 
ronment as suggested by Haas et al. (2004). 

Clearly, our data set is small and the results should 
finally be corroborated by larger samples from the en- 
tire Spitzer IRS data base. Nevertheless, since there are 
at most one (~15%) or two outliers in the distributions 
(Figures 2 & 3), the conclusions drawn from our 2x7 ob- 
jects may be generalized to 85% or at least 70% of the 
entire 3CR sources, i.e. the steep radio spectrum quasars 
and powerful galaxies with 10 26 5 W/Hz < Pi7smhz is 
10 29 - 5 W/Hz and 0.05 < z < 1.5. Also, the dispersion in 
the distributions is still large and needs to be explored fur- 
ther, possibly allowing for recognizing evolutionary trends 
among these complex sources. 
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